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A NEW GENERAL THEORY OF ERRORS. 

By William Edward Story. 

Presented April 13, 1904. Received May 5, 1904. 

I. Introduction. 

The ordinary theory of errors — the method of least squares — has 
been developed in three ways (methods of Gauss, Hagen, and Crofton) ; 
not only is each of these developments based upon assumptions to which 
not every one will assent, but at certain stages in each of the two later 
developments exactness is replaced by approximation through the neglect 
of certain quantities (or terms) regarded as very small (or productive of 
very little effect) in comparison with others retained, whereas this neglect 
is not fully justified, as it seems to us. 

The present paper is the outcome of an attempt to construct a theory 
of errors upon such simple principles as will be generally admitted to be 
necessary for the mathematical treatment of any such theory. The im- 
mediate object of the theory of errors is to determine the probability that 
the error made in the direct measurement of a physical magnitude shall 
lie between given limits. More complex applications of the theory 
(e. g. to functions of quantities directly observed and to quantities deter- 
mined by calculations based upon direct observations) can be connected 
with the solution of the fundamental problem just stated by means 
familiar to all who have used the ordinary theory. Of course we con- 
sider only accidental errors, assuming the observations to have been 
corrected for all systematic errors before being subjected to our investi- 
gation ; we suppose also that the observations were all made under the 
same conditions, so that they are equally reliable (their relative prob- 
abilities must be determined by the theory). 

II. Case op a Finite Number op Observations. 

Let 2 1; z 2 , z s , . . . , z„ be n observed values of a magnitude whose true 
value is Z, and let z denote an observed value in general. If A is any 
function of z and A lf A 2 , A s , . . . , A„ are the values of this function for 

*!> z 2> z z> ■ • • > 2»> respectively, we write 
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A 1 + A 2 + A s + ... + A n = [A], 

as is customary iu the theory of errors. Let z be the arithmetical mean 
of the observed values, so that 

(1) «o=i[*]. 

Let x be the error of an observation in general, x lt x 2 , x s , . . . , x n the 
errors of the n given observations, and x the arithmetical mean of these 
errors, so that 

{x = z — Z, a?! = z t — Z, x 2 = z 2 — Z, Xz = z z — Z, . . . , x n — z n — Z, 
1 l 

therefore, the arithmetical mean of the errors is the error of the arithmet- 
ical mean of the observations. Let £ be the residual of an observation 
in general, — that is, the excess of the observed value over the arithmet- 
ical mean of the given observed values, — and let &, £ 8 , £ 8 , . . . , £, be 
the residuals of the several observations, so that 

(3) £ = *-*„= (x+ Z)- (x + Z) =x-x ; 

therefore, the residual of an observed value is the residual of its error ; 
we may call it simply the residual of the observation. Let p k be the 
mean &-th power of the residual, so that 

(4) Pt = lffli 



in particular, 

(5) 






The positive square root of the mean square of the residual shall be 
called the mean residual (it is called the mean error in the ordinary 
theory) and shall be represented by /*,, so that 

(6) fi = Vj» 2 or p 2 = f*?. 

The ratio of the residual to the mean residual shall be called the relative 
residual and be represented by s, while the ratio of the mean /fc-th power 
of the residual to the £-th power of the mean residual shall be called the 
relative mean k-th power of the residual and be represented by p k , so that 



(7) 
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„ _ £ _ Pi 

* — ~> Pk — _ i> 



Po — — — h Pi — ~ — "> P2 = -j — 1- 
p n ix 



A little consideration will make it clear that no observations, however 
many, can determine the true value of the observed magnitude, nor, 
therefore, the errors of the several observations, — they can determine 
only what is equivalent to their arithmetical mean and the residuals of 
the observations. The fundamental problem of the theory of errors 
must, therefore, be modified, and may now be stated to be the determina- 
tion of the probability that the residual of an observation shall lie 
between certain limits according to the law of the long run, that is, in 
general, on the assumption that the distribution of errors over the whole 
range of possible error tends to follow a certain law, as the number of 
observations is increased without limit, and that the proportion of all the 
errors of any (finite) number of observations that fall between given 
limits is the same as if the number of observations was infinite. It is 
evident that the distribution of errors among different intervals in definite 
proportions implies an increase in the breadth of the intervals as the 
number of observations is diminished. The analytical theory of errors 
is based upon the supposition that the number of observations is infinite 
and that the interval for which the proportion of errors is to be deter- 
mined may be taken to be as narrow as we please. 

In connection with this statement o£ the fundamental problem it may 
be said that the determination of the value that may best be assigned to 
the observed magnitude is not explicitly a part of the fundamental prob- 
lem, although immediately connected with it and perhaps the most 
desirable result of it from the practical point of view, but the criterion 
for the best value is by no means evident. Perhaps it will be generally 
considered that the best value is the probable value, that is, the value that 
is just as likely to be exceeded as not in the long run, so that the prob- 
ability that the observed value shall fall short of the probable value is \ 
and the probability that the observed value shall exceed the probable 
value is |. Certainly the best value cannot be defined as the most prob- 
able value, that is, the value whose probability is a maximum, because 
it turns out that the probability is sometimes a maximum for more thau 
one value ; the probability may even have equal maxima for two or more 
values. 
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III. Case of an Infinite Number of Obseevations — 
Assumptions. 

The assumptions upon which we shall base the general theory of errors 
(upon which the ordinary theory has always been based in part, explicitly 
or implicitly) are these : 

I. Possible errors form a practically continuous sequence from a certain 
lower limit to a certain upper limit ; 

IT. The probability that the error of an observation lies between x and 
x + dx, where dx is infinitesimal, is <£(x) • dx, where <£(x) is an analyt- 
ical function of x developable by Taylor- s theorem throughout the whole 
range of possible error ; 

III. The probability that the error lies between given limits is independ- 
ent of the unit of measurement. 

Assumption I forms the foundation for II, and II implies that 
<jf> (x + h) is developable according to powers of h by Taylor's theorem 
for all values of x and x + h that lie within the range of possible error. 
The value of <f> (x) for any given value of x is called the relative fre- 
quency of the error x. 

The probability that the error of an observation lies between two given 
limits a (lower) and b (upper) is, evidently, 





(8) f+(x) 



dx. 



If x is the (algebraically) least possible error and x the greatest possible 
error, we have 

X 

(9) f<t> (x) ■ dx = 1. 

X 

If we express x in terms of the residual £ by (3), representing the lower 
and upper limits of possible residuals by | and |, respectively (which 
correspond to the extreme errors x and a;), and writing 

(10) 4,(*) = <Mf + *o) = ^(f), 
we have from (9), because x is constant, 

(ii) JV(3-#=i, 
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and the probability that the residual of an observation lies between 
a and b is 



(12) JV«)- 



d£. 



The value of i/r (£) for any given value £ may be called the relative fre- 
quency of the residual £, and \j/ (£) • d£ is evidently the probability that 
the residual lies between £ and £ + d£ ; that is, \f/ (£) • d£ is the propor- 
tion of all the residuals that lie between £ and £ + rff. In calculating 
the mean powers of the residual by (4), the power of each residual is to 
be multiplied by the number of times it occurs, and the sum of the 
products thus formed for all the residuals is to be divided by the number 
of observations ; if we assign all the residuals between £ and £ + d£ to 
the value £ itself (or replace them all by their mean), the multiplier 
corresponding to this value is proportional to \j/ (£) ■ d£, while the whole 
number of observations is proportional to (11). We have then 

f 

(13) Pt=f ***(£)•#, 

I 
and, in particular, by (5) and (6), 
I 

(14) j^(£)-d£= Pa =\, 
i 

1 

(15) ft*(.t)'dt=Pi = % 
i 

1 

(i6) fet(£)-d$= P2 = ,f. 

i 

Therefore, by (3), (14), and (15), 

(17) jx ■ 4(2) • dx = J(£+x ) ■ ,/,(£) • d£ = x 

2 i 

Because <]/ (£) is essentially positive for any possible residual £, formula 
(15) shows that the possible values of £ must include both positive and 
negative values, so that J < 0, < f, and is a possible residual in 
all cases. 
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IV. Effect of an Infinitesimal Source of Error. 

To determine the form of the function <//• (£) we consider the effect of 
the addition of an infinitesimal source of error to the original source 
(however complex that may be), that is, of a source of error that pro- 
duces infinitesimal errors with great frequency and finite errors of even 
moderate magnitude with infinitesimal frequency. We employ the same 
notation for all functions and quantities pertaining to the action of the 
infinitesimal source alone as for those pertaining to the original source 
alone, distinguishing them by an accent, and denote changes due to the 
addition of the infinitesimal source by the symbols for the functions and 
quantities changed with a 8 prefixed. 

A given error x produced by the simultaneous action of both sources 
is the result of a certain error x' produced by the infinitesimal source and 
an error x — x' produced by the original source. The probability that 
the error produced by the combination shall lie between x and x + dx 
is, then, 

x' 

|> (a:) + H (x)] dx = f<t>(x- x>) • <f,' (*') • dx ■ dx', 
from which follows 
(18) <t> (x) + 8<j> (x)=J<t>(x- x>) • 4>' (*') ' dx'. 

x> 

For the arithmetical mean of the error produced by the combination 
we have 

~x-\-5x Ic 1 

x + 8*0 = fx- dx- f<t>(x - x') • 4>> (x<) • dx' 

x+Sx x] 

X x< 

= fx-dx- C<j>(x — x')-<f>'(x')-dx' + x-cj>(x)-8x-x-<j>(x)-8x. 

x x[ 

It is to be observed that, in whatever manner these integrals are 
evaluated, x — x' must be taken over the whole range of values from x 
to x. The limits of possible error (x and x) cannot be fixed a priori ; 
some may think that large errors or residuals are impossible by the nature 
of accurate observation, but we prefer to regard the limits as determined 
by the frequency function <j> (x), namely, we assume that the value of 
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<f> (x) (which is positive for any possible error x) is at each limit ; * 
if <f> (x) is not for any finite value of x, all values of the error are to be 
regarded as possible. The results will show that <j> (x) = for x = ± oo, . 
and, in general, <f> (x) has a very small value for any large value of x. 
We have, then, 

(19) <f> (x) = and <f> (x) - 0, 

so that the last two terms of x + 8x vanish, and, on reversing the order 
of integration, by (17) and (9), 

x' x+x' 

x + Sa? = / <f>' (x 1 ) ■ dx' ■ I x • <j> (x — x') • dx 

xi x+x 1 



X' z 

= j V (x>) ■ dx' ■ C(x + x') • <j> (x) • dx 

X[ X 

x' 

r= J (x + x') <j>' (x>) dx' = x + x ' ; 



so that 

(20) Bx = x >. 

By (10), as applied to both component sources and their combination, 
we have 

^(x') = ^^ + x ')=^^), 

< i> (x-x')=4>(£-£' + x s> - oV) = <p(£-e- *„'), 

<f> (x) + 84, (x) = (4> + o<£) (( + *„) = (^ + fy) ($ - x ') = 

4, a - xj) + ty((- x '), 

* The relative frequency ^ (|) of a possible residual £ must be positive ; if | 
varies continuously from zero, either through positive or through negative values, 
the limit | or j must be reached before ty (|) becomes negative ; but if \f/ ( £) passes 
through zero without changing sign as { passes through a certain value, — a value 
that corresponds to a point of contact of the curve represented by y = t\i (|) with 
the axis of abscissae, — that value of £ is not necessarily a limit to the possible 
values ; such a value of | could never occur as a residual, although values on both 
sides of it might occur. However, if the conditions under which the observations 
are made change continuously in such a manner that the curve y = $ (£) assumes a 
form in which it touches the axis of abscissae, this must be the limiting (extreme) 
form, provided a continuation of the change would produce a curve having negative 
ordtnates in the neighborhood of the point of contact ; namely, on account of the 
permanence of the relation (14), the limits J and J cannot change discontinuously. 
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for any given value of £, because the arithmetical mean of the error cor- 
responding to the function <j> + S<j> is x + Sx = x -f- x '. Therefore, 
(18) becomes 

f 

or, when f is replaced by f + a? ', 

r 

(21 ) $ ({) + ty (0 =f*Q- £') ■ V (£') ■ d£' 

'(' 

for any possible value of £. Now, by assumption II, <f> (x) is develop- 
able by Taylor's theorem and, therefore, so is if/ ($), that is 

o 
therefore, by (13), (14), and (15), 

* (o + ^ (i) =2 ft (- *)* ;n ^|!P /> f (O • #' 

i' 

- ^ ; * ! 3i* 
o 



-♦(e+^-^ft-s 



so 



that 



(22) *<0=Jt(-D'?i^ 

2 



Also, by (13) and (21), 
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Reversing the order of integration, observing that £ — £' must take 
all values from | to | and that, by (19), >f> (S ) = and \p (|) = 0, we 
have 

f T+(' 

p m + *p m = JV «0 ' <% 'ft m ■ t (t ~ O • « 
£' i+f 

r f 

=J ,// (o • «*F • J*(£ + £T • -Ml) • « 
£' i 

= JL ( " ) J r ' * ( *> ' *' '/*" * • ^ ^ ' * 
o i' .« 

2 

f Tfl\ 

where ( , ) denotes the coefficient of x* in the development of (1 + a;)™ 
by the binomial theorem, that is 

Therefore, 

m 
(24) Bp m =^> (™\ Pk'p m -k for 2 < m. 

2 

In particular, for the smaller values of m, 

Sp = 0, Spi = (because p = 1 and ^ = 0), 

8?> 2 = £> 2 ', 5j» 3 = p 3 ', 8p 4 = 6 p^pj + pi, 

lp h = WpsPi + 10p 2 p 3 ' + p 5 ', 

Sp 6 = 15p i p 2 ' + 20p 3 p d ' + IbpzpJ + p e ', etc., 

and, therefore, 

Pt = ty* > Pz = 8p« > P* = ¥4 - 6 />2 • 8p 2 , 
K = fy> 5 — 10^> 8 • Sp 2 — 10 p 2 • Sp s , 

Pe = $P« ~ ^ 5 Pi ' ¥2 - 20p 3 • Sp 3 — 15 p 2 (o> 4 — 6j» 2 - 8p a ) 
= 8p 6 — 15p 2 • &p t — 20p s ■ Sp s — 15 (p t — 6p 2 2 ) 8p 2 , etc. 
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Because the right member of (24) is linear in the p' 's and pj has the 
coefficient 1 in the expression for hp m but does not occur in the expres- 
sion of any Sp with suffix less than m, p m ' can be expressed as a linear 
function of the 8/?'s with suffixes not greater than m whose coefficients 
are polynomials in the p's. The expression for p 6 < in terms of Sp 2 , 
8p s , Spi, and Sp 6 shows that the coefficients of the Sjo's in the expres- 
sions of the p' 's are not generally linear. If we substitute the expressions 
of the p' 's as linear functions of the Sp's in (22) we obtain Si^ (f) as a 
linear function of the Sp's whose coefficients are polynomials in the p's 
multiplied by derivatives of <]/(£) with respect to £ ; therefore, i/' (f) is a 
function of £ and the p's with suffixes as great as 2 alone, say * 

(25) $ (() = F (£, p s , p 3 , PtoPs,. . .)• 

This result confirms the statement made above that all theory can 
determine is the probability that the residual shall lie between given 
limits. 

V. Change of the Unit of Measurement. 

If we reduce the unit of measurement in the ratio 1 : A, the numerical 
values of all measurements are multiplied by A, so that, for given obser- 
vations, x, x , i, d$, p k (for 2 < k) have to be replaced by A x, A x , A £, 
A • d£, X k p k , and, by assumption III, the probability that the residual 
lies between A $ and A £ + \ • d£ with the new unit is the same as that 
it lies between £ and £ + d£ with the old unit ; therefore, 

F{\£, A>, X s p s , A*/>o . . .) • \d£ = F(£, Pz,p s , Pi ,...)- d£, 
that is 

(26) F(X£, A> 2 , A>, k*pi, . . .) = I F(£, p„ p s , p t , . . .) 

for all values off, p%, p s , p it . . . The relations (25) and (26) express 

that>i^(f) is a homogeneous function of £, \Zp%, Vp s , Vpa etc. of de- 
gree — 1, that is, with the notation of (6) and (7), 

C 27 ) <A (£> = -/(«> Ps, P*> Pi, ■■■)■> 

r 

for convenience, we shall sometimes abbreviate \fr (£) to ij/ and 

/(*, p„ Pi , p 6 , . . .) to f(s) or/. 

* This conelusisn implies that every infinitesimal change in tji (|), for a given 
value of |, can be effected by the addition of an infinitesimal source of error, — 
which ought, perhaps, to be stated as a fourth assumption lying at the base of the 
general theory. 
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VI. Differential Equation of /(s). 

Formulae (22) and (24) give the variations of tp (£) a,ni p m due to the 
addition of any infinitesimal source of error ; the mean square pj of the 
residual due to such a source is an essentially positive infinitesimal ; in 
determining the mean powers of the residual for such a source by (13) 
we may confine our attention to values of the residual that lie within 
very narrow limits which, in accordance with the remark made in con- 
nection with (15), include 0, so that we need consider only very small 
values of $; the terms of the integral (13) for greater values of k are 
infinitesimals of a higher order than the corresponding terms of the 
integral for smaller values of h, so that the integral p t ' for any value of 
k greater than 2 is an infinitesimal of a higher order than pj. We 
specialize the infinitesimal source of error by assuming p k ' = for 
3 ^ k* We have then, by (24), 



fym = ( 2 ) P m ~ 2 ' P * i0T 2 < m ' 



from which, together with (6), follows 

p./ — &p i = 2/x- S/x, 
so that hp m = m (m — 1) • p m _ 2 /i • Sii, 

while (22) gives 

H (c) = \n g^ • = g^ •/*•*>/*• 

We denote partial differentiation by the round 3, regarding i/< as a 
function of £, p 2 , p s , Pi, ■ ■ ■ and/ as a functiou of s, p 3 , p 4 , p 5 , . . . , so 
that, by (7) and (27), 

9l_}_9f 3v__isy 

9i ~ /x 2 3s 3f 2 ~ ^ 3s 2 ' 

The variations denoted by 8 were all taken for a definite value of £; 
therefore, by (7), (27), and the formulae just written, 

S/U. 



= 8(-)= — -j-8/* = -« 

ty = 3p „ 8 M = ^ -j, 



f* 



* This assumption does not in any way affect the generality of the results 
obtained, because p k ' enters into <|< (|) only through the variations of the unac- 
cented p's. 

VOL. XL. — 12 
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S/ = 80. ^r) = p • ty + + • 8/* = \% + f\ j> 
Sp* = A (# — 1) • Pt_2 p-* -1 - 8/i, 

Sp* = s(^)=%-*-^-Sp. = [*(*-l)- P *_ 2 -*- Pi ] 
therefore, 






= ^8s + 



2*£w' 8/,t 



5« -""^i cW 

3 

={-g+£[»0-'>~-w£}£- 

Equating the two expressions for 8/ thus found, omitting the common 

factor — > we obtain 
P- 

< 28 > ^+'S +/ =?* C * ( *- 1 ^-*^^' 

a differential equation for/ which may be written 

< 29 ) a(l + ^)- 8 -^ 

where 

(30) 7" = i]^ [*(*-!) P*- 2 - *P*] I; 

is a differential operator that aetfs wjoo?! f on/y through the p's with suffixes 
as great as 3 ; the first part of that term of T for which k = 3 vauishes 
because p t = 0. We have omitted the upper limit to the sum in T be- 
cause that depends upon the number of p's involved in f, a number as 
yet unknown. 

VII. The Operator T. The P's and S's. 

The operator T is so important for the theory that it deserves par- 
ticular study. By a constant we shall henceforth mean a quantity that 
is independent of the p's with suffixes as great as 3 as well as of s. If 
u and v are any two operands and O is any constant, we have 



(31) 
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T (u + v) = Tu + Tv, 

T- Ou= G- Tu, 

T- uv = u- Tv + v ■ Tu. 



If u is a function of the p's with suffixes not greater than m, Tu is also 
a function of the p's with suffixes not greater than m ; if u is a poly- 
nomial in the p's, Tu is a polynomial in the p's of the same degree as u ; 
in particular, if u is a linear function of the p's with suffixes not greater 
than m, Tu is a linear function of the p's with suffixes not greater than 
m. If 3 < m, 

(32) T Pm = \ m (m — 1) • Pm _ 2 — \ m • p m . 

Therefore, if m is a linear function of p m and the p's with suffixes less 
than m, where 3 < m, T" u is a linear function of p m and the p's with 
suffixes less than m in which the coefficient of p m is the coefficient of p m 

(m\ n 
-- , where n is any positive integer or ; more 

generally, if G(T) denotes any polynomial in T (linear function of 
powers of T with constant coefficients), G (T) • u is a linear function 
of the p's with suffixes not greater than m in which the coefficient of 

p„ is the coefficient of p m in u multiplied by G I — — ) ; in particular, 

G ( T) ■ u will not involve p m if G ( T) contains T + — as a factor, and 
only then. Formula (32) holds also for m = 0, 1, and 2, giving 

Tpo = • po = 0, T'pj = — J Pl = 0, yp 2 = p — p 2 = 0, 
because p = p 2 = 1 and p x = 0, by (7) • therefore, 

(7*+ f) po = 0, (T + i) Pl = 0, (T+ f) p 2 = p 2 = 1 ; 
so that G(T)- p m does not involve p m if ©(T 7 ) contains the factor 

Tit 

T + — , where m is or any positive integer excepting 2, and G(T) • p 2 

does not involve p 2 if 6r ( T) contains fasa factor. Of course, poly- 
nomials in T are commutative with each other, so that it is of no conse- 
quence what position in a product G(T) the annihilating factor occupies. 

It is convenient to employ the usual symbol.y" 1 to denote the product 
of n factors diminishing by successive units from y, whatever meaning y 
may have (it may be a symbol of quantity or of operation), so that 

(33) ^ =y (y - 1) (y - 2) (y - 3) . . . (y - n + 1) 
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if n is any positive integer ; in consequence of this definition, 
(34) 



/y» = y, 



,(») _ 



r 



y — n 

for 1 < 'n ; we take the last formula as the basis of the definition of y'" 1 
f or n = 0, — 1, — 2, — 3, . . . , successively, and thus obtain 

1 



(35) y<°> = 1, 



,,(-»> 



(y + 1) (y + 2) (y + 3) . . . (y + n) 



for any positive integer n (that is, for any negative integer — n). With 
this notation we have 

y> = yW • (y — 4 - )<"- j » = (— 1)" • (— y + n — 1)'">, 

(y + » - i) w - (- i) n ' (- y) (n, > 

w ! 



(36) 



(n - i) : 



lCK=(VW;::> 



where each of the letters n and i represents any positive integer or ; 
also, 

r(20! = *"»!(»-i) w 
( ' \ (2 * + 1) ! = 2 2i+1 • t ! (» + |)<' +1 », 

because, evidently, 2' i ! is the product of the even numbers from 2 to 2 i, 
2' (i — -iy > is the product of the odd numbers from 1 to 2 * — 1, and 
2' +1 (i + i) (m) is the product of the odd numbers from 1 to 2 i + 1. 

Furthermore, we have 

■ 
(y _ 1)<») — (y — I)*— 1 ' • (y — n )=y in) — n ■ (y— 1)<"- 1 > 

. w (3) . y»-3) + . . . 



— y(») _ w . yn-D + n ' 



.(»> • „C-2! 



(38) 






«<*> • y'"-*), 



where n is any positive integer or 0. 

From (32) it is evident that any polynomial in T changes p with any 
even suffix into a linear function of the p's with even suffixes and p with 
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any odd suffix into a linear function of the p's with odd suffixes. There- 
fore, if u is a linear function of p 2m and the p's with even suffixes less 
than 2 m, where 2 < m, (T + m) u does not involve p 2 „, (T + m) ,2) it 
involves neither p 2m nor p 2m _ 2 , (^"+ m ) m u involves neither p 2m , p im _ 2 , 
nor p 2m _ 4 , and so on, and (T + m) (m_1) M is a constant (or linear function 
of p and p 2 ) ; therefore, by (36) 

( (T+ mY m+1) u = (T+ mY m - 1) (T+ 1) (2) u 

(39) \ 

that is, (T 7 -)- m) (,n+1) annihilates every linear function of the p's with; 
even suffixes not greater than 2 m, where m is any positive integer or • 
(this is true for m = or 1 because the operator contains T as a factor). 
But ( T + m) |m+1) annihilates no linear function that involves a p with- 
an even suffix greater than 2 w or an odd suffix as great as 3 ; namely^ 
if p k is the p with greatest even suffix or greatest odd suffix (other 
than 1) that is involved in the linear function u (with a coefficient dif- 
ferent from 0), the coefficient of p h in (T + m) im+l) u is the coefficient 

of p k in u multiplied by I m — - ] , whose factors are all negative if 

k is an even number greater than 2 m and all odd multiples of f if h 
is odd. 

Similarly, if u is a linear function of p 2m +i a "d the p's with odd suffixes 
less than 2 m + 1, where < m, 

(40) ( 7 1 + m + i)"" +I » » = 0, 

so that (7"+ m + £) (m+1) annihilates any linear function of the p's with 
odd suffixes not greater than 2 m + 1 ; but (T + m + %) lm+l) annihilates 
no linear function that involves a p with an odd suffix greater than 
2 m + I or an even suffix. 

Applying (39) and (40) to p 2ra and p 2m+1 , respectively, we have 

( T + ™) ,m+1) p 2m = 0, (T + m + i) ,m+ » p 2m+1 = 0, 
so that, by (34), 

T(T+m-l)™ P2m = -m- (T + m - l) im) P2m , 
T(T+m- ir ] P2m+ i = -(m + i)-(T+m- *)*"» p 2m+1 , 



(41) TQ 2:n = -m- Q 2m> T Q 2m+1 = - (m + *■) ■ <? 2m+ i, 
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if we write 

(42) (T+m- 1)"»» P2 „ = Q 2m , (T+ m - i)™ P2m+1 = Q 2m+l ; 

and, evidently, formulae (41) hold for any constant multiples of Q 2m and 
Qzm+i, respectively. 
Now (32) shows that 

(T + m) • p 2m = 2 m (m - $) • P2ril _ 2 

and, by continued application of the same formula, 

(43)* (T + my» P2m = 2* m<*> (m - £)<*> • p 2m _ a 

for < h < m — 1 ; in particular, 

(T+ m) , "- I V 2 ,„ = 2"- 1 m' m -» (m — f) <ra -»p 2 
and 
(44) ( T + «)<"" p 2m - 2"* • m ! (m - £)"»> Po , 

because, by (34) and (36), 

m ! = mf m) = mf m - 1} • 1, 
(m - £>""> = (m - i) 1 " 1 -" • i, 
(T+ \)p 2 = p 2 = p ; 

so that (43) holds also for k - m. Then, by (42), (38), (43), and (36), 



(45) 



o 

m 

= 2- (m - *)» • » !^(- 1)' (j J F 



Pv 



■ (•-*)• 



for < m. 

Similarly, we have, from (32), 

(2 1 + ot + i) ■ p 2m+ i = 2m(m + i) • pt^! 
and, by continued application, 

(46) (T+m + *)'*> • P2m+1 = 2* «« (m + £)<*> • p^^., 

for ^ i ^ m. Therefore, by (42), (38), (46), and (36), 

(47) 



= *•« (» + i)«-H> • » ! j? (- i)<(™) 2^^ 

for <* m. 
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If we write 
r 



(48) 



m 



?w_ i/ l) VJ 2W 



PSi+1 



!' +1 (* + i) (i+1> ' 
we have, by (45) and (47), 

. / &» = 2" (« - i) (m) « ! -P.. and 

I <? 2m+1 = 2- +1 (m + *r«' » ! iWi. 
and, because P 2m and P 2M+ , are constant multiples of Q 2m and ^ 3m+1 , 
by (41), 

(50) TP 2m = -m- P 2m , TP 2m+1 = - (in + *) • P 2 » +1 , 

for all values of m. For the P's with the smaller suffixes we have, from 
(48) and (7), 

(51) f Po=1 » ^1=^ = 0, P 8 = -Jp 3 , P 4 = -1 + Jp 4 , 
\P 5 = -|p 3 + APs, Pe = -2 + p*-T 1 ]fP« etc.; 

so that equations (50) give 

(52) rP = 0-P„ = 0, 7'P 1 = -£P 1 = 0, TP 2 = -P 2 = 0. 

Because P^ is a linear function of the p's with even suffixes up to 2 m 
in which the coefficient of p 2m is different from 0, and P 2m+ i is a linear 
function of the p's with odd suffixes up to 2 m + 1 in which the coefficient 
of p 2m+1 is different from 0, p 2m can be expressed as a linear function of 
the P's with even suffixes up to 2 m and p 2m+1 as a linear function of the 
P's with odd suffixes up to 2 m + 1. In fact, if we multiply P 2m or P 2 m+i 
by (_ i) m ( ) and sum for all values of m from to h, we obtain 

£<-K:)^4 ( - iy *^£ ( -' r (X) 

Pat 



2*(i -*)<*>' 

* m 



^-^Ciw^-w^iX^oc 



P 2 lt+1 



2*+' (k + i)'*"™' 
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because 

k 



£.<-<) (7) =(?)§.<- >Ki-') 



therefore 



= (~ 1)' 

= if i < k and 
= (-1)* if » = *; 



C)i<-<:0 



(53) 



P, = 2* (1 - i) ! *»2 m (- 1 )'" (I) P *» for ° < *» 

p a+1 = 2W (i + i)* +1, ]§ m (- l) m ( *)-P-+i for < *• 



o 
As a consequence of (50) we have 

(54) ^(TVi^^-^-P™, 

if G(T) is any polynomial in T. The ease with which the result of 
applying any polynomial in J" to the P's can be obtained by means of this 
formula makes the introduction of the P's in place of the p's very advan- 
tageous. We shall henceforth assume that all functions of the p's have 
been expressed as functions of the P's by (53), so that the degree of any 
polynomial in the -P's with suffixes as great as 3 is the same as that of the 
corresponding polynomial in the p's. Equations (50) show also that the 
operator T, defined by (30) in terms of the p's, can be expressed in terms 
of the P'srthus : 

'(55) T^-i^kP,^. 

We assign to each P a weight equal to its suffix (leaving out of account 
P x and P 2 , which are and, therefore, nowhere occur) and define the 
weight of any product of P's to be the sum of the weights of its factors ; 
any constant may be regarded as a multiple of P = 1, and, therefore, as 
of the weight 0. The weight of any such product as P 3 n s P 4 TC i P 6 »5 ... is 

(56) 3% + 4»4+5kj + ... = w, 
say, so that, by (55) and the last formula of (31), 

iff) 

(57) T • P s n * P/ 4 P 5 " s . . . = - ^ • P 3 % P 4 " 4 P"" . . . 
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and, evidently 

(58) G (T) • P 3 n » P 4 "< P 5 " 6 ...= g(- I V P s % P 4 re * P^ ... 

where G(T) is any polynomial in T. More generally, if II is any 
isobaric polynomial of weight w in the P's, that is, any polynomial whose 
terms are all of one weight w, 

(59) «?(T)-n=ff(-|)-n. 

Therefore, if the coefficients of II are not all 0, (? (y) • II is or is not 

w 
according as G {T) does or does not contain T+ — as a factor. In 

particular, (T + tn) {m+1) annihilates every polynomial in the P's whose 
terms are all of even weights not greater than 2 m and annihilates no 
polynomial of other description, while (T+ m + ^y m + 1 > annihilates every 
polynomial whose terms are all of odd weights not greater than 2 m + 1 
and no others. By means of formula (59) any polynomial in T can be 
applied to any polynomial in the P's ; namely, any polynomial in the 
P's can be expressed as the sum of isobaric polynomials of different 
weights and the polynomial in T can be applied to each isobaric part 
by (59). 

The P's can be determined from the function /(s) ; namely, from (7), 
(13), and (27) we have 

(60) Pk = g = P~ k + (0 ■ di = fs* ■ f(s) ■ ds for ^ k, 

i 1 

and obtain, by the substitution of these values of the p's in (48), 
s ? 

(61) P 2m = Js 2m ■/(.) ■ ds, P 2m+1 = JS^ X -f(s) . ds, 



where 



(62) 



.«2'+l 



^»+i -£jf~ ^\i)W(i- + iV" ! '+i> ; 



2 !+1 (T+ iy 



for < m. For the S's, with the smaller suffixes we have 
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(63) 



'S = l, S t = s, S 2 = 1 — s 2 , S 3 — s — £ s g , 
# 4 = 1 - 2 * 2 + ^*, -S. = «-§j , + A'*. 
[S„ = 1 - 3 s 2 + s 4 - T 'j s 6 , etc. 



VIII. The Form of /(«). The i?'s. 
To simplify the integration of equation (29), we put 
/(*) = (s) ■ u, 
where u is an integral of the equation 

say 

u — e %' 
where e represents the base of natural logarithms, so that 

(64) f(s) = 0(s) -e-f. 

Writing $ as an abbreviation for (*), we have 

rf = L^ 2 '^~ ( } _T*' 

so that (29) becomes 

T3 2 3(91 «' „,„ « 2 

_* ! 
because e 2 is independent of the p's (and P s) ; therefore, 

<«) S — B— r * 

which is a differential equation for (s).' 

The assumption II that <£ (a;) is developable by Taylor's theorem 
throughout the range of possible values of x implies that f (f ), f(s), 
and (s), are similarly developable ; we may therefore put 

(66) Hs)=S,%s\ 



where A { (for each value of i) is a function of the P's alone. This de- 
velopment implies that we shall approximate more and more nearly to 
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the true form of 6 (s), even if its complete expression is an infinite series, 
the more terms of (66) we take into account. Taking any finite num- 
ber of terms, say for values of i up to k, substituting in (64) and the 
resulting expression for/(s) in (61), we obtain 

(67) P "=%jtJ siR <"' ds > 

s 

where 

(68) R m = S m - eT% for ^ m. 

We see from (67) that, to the degree of approximation that corresponds 
to k + 1 terms of (66), the P's are linear functions of the k + 1 A's, 
whose coefficients are definite functions of the limits s and * alone, because 
the (S's and, therefore, the P's are independent of the P's. Therefore, 
not more than k + 1 of the P's (including P , P x , and P 2 ) are linearly 
independent ; the solution of such equations (67) as correspond to linearly 
independent P'a gives the k + 1 A's as linear functions of these P's, 
whose coefficients are definite functions of * and * derived from the inte- 
grals involved in (67) ; these expressions of the A's will involve a certain 
number of arbitrary constants if the number of linearly independent P's 
is less than k + 1, but that is of no consequence. If the P's are sub- 
jected to infinitesimal variations, — say 8P ra is the variation of P m , — 
the A's, f(s), and (s) will undergo variations, which may be repre- 
sented by SA e (for A,), 8/(s), and 80 (s), respectively. Then, by (66) 
and (64)", 

^-A Si S 2 

80 W =\ if $i > 8/(S) = W {S) ' ^ 

and, by (61) and (68), because f(s) = for each of the limits s and s, 

s" s" s 

8P,„ =J-S m ■ 8/(.) • ds =JR m ■ S0 (») • ds =2 ti -jff? R m ■ ds ; 

that is, the relations between the 8P's and the 8^4's are the same as if 
the coefficients of the A's in (67) were constants (whether £ and s are 
functions of the P's or not) ; therefore, in determining the A's as func- 
tions of the P's that satisfy the differential equation (65) we may regard 
them as linear functions of the P's with constant coefficients. The 
necessary and sufficient conditions for these linear functions are those 
implied in (65) together with 
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(69) A, = for K < i 

(to the degree of approximation to which we restrict oar calculations). 
Substituting the expression (66) for 6 (s) in (65) we obtain 



<J* (i - 2) ! 4(»-l)! A»! ' 



that is 



" 1— 

£0- 



;^,-2-rj,Jij = o 



for every possible value of s, where the first term of the multiplier of — 

is to be omitted for i = k and i '. = k — 1 and the second term goes out 
of itself for i = 0. Therefore, the necessary and sufficient conditions 
for the satisfaction of (65) by (66) are 

A i+i = 2 (t+ lV A, for ^ t ^ k - 2 
and 

( T+ I J • ^, = for » = k or K — 1 ; 

that is 

(70) J '" 2 ( r+ T)' 4! for 2 <*< " 

and 

(7i) (r+i)-^=o 

for -4, with the greatest even suffix and A, with the greatest odd suffix 
that occur in (66). On replacing i in (70) by 2 i and 2 % + 1, respec- 
tively, we have, for the A's with even and odd suffixes separately, 

A v = 2 {T+ i - 1) • J 2 ,_ 2 for l^i<l> 

A+i = 2(7"+ *-!)• A_i for 1<*"<^, 
and, by repeated application of these formulae, 

(72) A 2 , =2'(T+i- 1)« -J, for l^i'^ 
and 

(73) J w = 2<(7*+ ■-*)«• .4 for 1 = ^^-1; 
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that is, the A's with even suffixes are all derived from A and all the A's 
with odd suffixes from A 1 by the application of definite polynomials in T. 
We may, then, determine the forms of the A'a with even suffixes and of 
those with odd suffixes separately ; that is, we may suppose our calcula- 
tions restricted to a degree of approximation that corresponds to a devel- 
opment of 6(s) in which occur even powers of * up to s 2 <* and odd powers 
up to s 2 ^+ 1 , without assuming any particular relation between a and /3. 
Then (72) holds for l^i^a, (73) holds for 1 ^ i ^ /?, while (71) is 
expressed by 

(T+ a) • A ia = 0, (T+ p + |) • A 2li+l = 0, 

that is, when A^ a and -42,3+1 have been determined by (72) and (73), 

respectively, 

(74) (T + a)(*+!) • A = 0, (T+ /3 + £)P+i> • A x = 0. 

The only conditions imposed upon A and A x so far are, then, that they 
shall be linear functions of the P's with constant coefficients that satisfy 
the equations (74). The only linear functions of the P's that are of 
a given weight w are constant multiples of P w , while, by (54), 

(7'+a)(«+D- J P„=^a--j -P., 

(T+p:+0?+V-P m =((3 —J -P„; 

therefore, in order that A (as a linear function of the P's) shall satisfy 
the first equation (74) it is necessary and sufficient that the suffix w of 
any P involved in it shall satisfy the equation 

\(«+l) 

= 0, 






that is that w shall be an even integer not greater than 2 a ; and, in order 
that A x (as a linear function of the P's) shall satisfy the second equation 
(74) it is necessary and sufficient that the suffix w of any P involved in 
it shall satisfy the equation 

that is that w shall be an odd integer not greater than 2/3+1. We 
may, then, write 



(75) A Q =^ k G u P u ,A 1 =2 4 

o 
where the C"s are constants. 






(76) 
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On applying (72) and (73) to the first and second of equations (75), 
respectively, we obtaiD, by (54) and (36), 

- a a 

A *=% k °» • p * 2i (»-*- i ) (i '=2* cs * /v (-!)•' 2<km > 



*■ 

Substituting these values of A 2i and -4 2 ;+i in (66) we have 

o o 

+2^ k °^ p ^ '2iS~ ^ (2.- + i)i ^ 



or, replacing (2 1) ! and (2 t + 1) I by their values in (37), observing the 
last line of (36), and employing the notation of (62), 



■^jfc+l ^!itl "2*+l ) 



(77) 6 («) =^ k G 2i P u S» +^ k C 2k+1 P 2i 

o o 

again, by (64) and (68), 

a £ 

(78) f (*) =z j ^^ l . u " 2 * 2 * ~^~ ^Jk k+1 i * +l 2 *+ 1 * 

o o 

Substituting this expression of/(s) in (61) we obtain 

Plm— ^,.C'tkPtk' I "2mRu'd S + ^^ @1k+\P 2*+l' I "2m^2*+l" s > 



(79) 



s 

9 



^ > !«+l = ^ i ^»°» / "2m+l-"2*'^+^V Cz*+l-«2*+l" l"4»i+l"2t+l'<"' 



IX. The Coefficients of /(s). Conditions fob the Limits. 

In all that contributed to the determination of the form of f(s) as 
given by (78) no assumption has been made with regard to the values of 
the P's with even suffixes not greater than 2 a and the P's with odd 
suffixes not greater than 2/3+1 [excepting that, by (51), P = 1, 



(80) 
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Pi — P z = 0], — in fact, they have been assumed to be entirely arbi- 
trary ; therefore, the first of equations (79) for < m < a and the 
second of these equations for ^ m < /? must be satisfied independently 
of the values of the P's involved ; that is, for any one of the values of m 
in question, the coefficient of any P in the right member of the corre- 
sponding equation must be excepting the coefficient of the P that 
occurs in the left member, which must be 1. Therefore, 

Gim I S 2m Rim ds =1 for < m < a, excepting m = 1 ; 

_s 

s 

CL+i / Sam+i fi im+ i ■ ds = 1 for 1 < m < /3 ; 

£ 
J 

j S im R n • ds = for 0<m<a, 0<£;5 a, k ^ m; 

s 
I S 2m R u+1 • ds — for < w <j a, <£;?/? ; 

(si) if_ 

/ *^2m+i -#2* ' rfs = for < m < /2, ;? £ < a; 

s 
s 

/ Sim+i ft^+i • ds = for < m ^ /3, < ^ < /8, A ^ m. 

From the first, second, third, and fourth sets of equations (81) we ought 
to exclude cases for which k = 1, 0, 1, and 0, respectively, because 
Pi = Pi = and the terms corresponding to these cases do not occur in 
(79) ; but, by (68), 



(82) 



S S s 

js m R„ ■ ds =Cs m S ■ e~* ■ds = Cs k R m - ds, 



so that the cases that are to be excluded from any one set of equations 
(81) will occur (with an interchange of m and k) in other sets, excepting 
only the case of the integral 
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s s 

/ S% Ri • ds = j Si R 2 



ds, 



which does not occur as the left member of any equation. The four sets 
of equations (81) can be represented, by virtue of (82), as the one set 



s 

(810 /. 



SLB„-dt = 



for any two different numbers m and h from among the even numbers 
not greater than 2 a (including 0) and the odd numbers not greater than 
2/? + 1, excepting m = 1, k = 2, and m = 2, £ = 1. The equations 
(80) determine all the constant coefficients G that actually occur in (78) 
(Ci and C 2 do not occur), while equations (81) express relations that 
must exist between the limits s and s for any observations that conform 
to the theory. 

To determine the values of the integrals that occur in (80) and (81), 
we consider the functions S and R defined by (62) and (68). We have 

9s ~ J*C ' \ i J 2'- 1 (t - *)«> 



m 



i) 



(i) 



m—1 



m 



that is, 

(83) ^ = -2«'Vi, ^±i=S s 

for < m. We have also, by (36), 

m— 1 m 

^ = 2/~ 1} \ •■■ ) 2'(.-l*)«» + 2i (_ 1} v - 1 J^t-ir 

n-l 



= *~.-2 ( <-vC"7> 



c 2i+2 



+1 (»■ + !)«+« 
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(84) = S 2m _ 2 — s • Sin-i for 0^»i, 

m 

s ■ *. =^(- iy (7) [2 m + 1 - 2 (» - ,-j] 2i+1 ( ; 2 ; +1 i)(i+1) 
= (»»+!)• ^ - S-2J/- d'(" T % +1( ;; +1 i)(W ' 

= (2 ro + 1) • £ 2m+1 -2m- .S^-i , 
that is, 

2 m 1 

(85) fl^ = £__.,S^ +___.,. ^ forO^m. 

From (68), (84), and (85) we obtain 

{Rim = Rlm-1 — * " -Ram-l for < TO, 
2m 1 

i?2m+1 = 2^TT ' *«-» + 2^T1 " ' ' *- for ° < m ' 

while from (68), (83), (84), and (85) follow 

3R« s 2 s 2 

-gf = -(2»- Sam-i + « ■ £,„,) e J = -(2«+l) S 2m+ i • e~* 

(87) =-(2 ffl +l)'4 +1) 

— ^f^ = (S 2m — s " ^a m +i) e * = ^ 2m+2 • e~ 5 

(88) = i? 2m+2 . 

From (83), (87), and (88) follow, after simple reductions by (84)-(86), 
9 

■q s i m ' ^an-i Rik—k' S 2m i? 2 *-l} = (»» — *)• /Sj, m i? 2 * , 

3 

^~ \»J2m+l -"24 — "2m -»2*+l} = 2 (TO — K) ' iSjm+l -""2*+l > 

5 

^|2»»' »%,»-! ^ 2i+1 + S 2m R 2k } = (2 m - 2 £ - 1) • S 2m B u+1 . 

Therefore, letting [ ] denote the excess of the value of the enclosed 
function of « for s = s over its value for s = s, we have 



(89) / 



m 



S 2m R 2k ■ ds — m _ k [^2«-l #2*] £ — m _ k L S 2mRu-l]l 
VOL. XL. — 18 
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for k ^ m, 

(90) 

s 
/ ^2m+i -Rsw-i ' * = 2 (m _ jfc) [^ 2 ™+i ^2*] £ ~~ 2(m — k) ^ 2m ^ 2S+1 ]s 

for k ^ m, and 
(91) 

J ^^+1 " * = 2m-2k-l t^-i^X + 2»^2^ri [^aft 

Also, from (83), (87), and (88) follow 
9 

W~ (<J27« ■Kim—l) = "2m -"-2m — 2m* 02m-X -"2m— 1 • 

3 

p~ (»$Jra+l #2m) = — (2 *M + 1) • S 2m+i i? 2m+l + ^2m Am • 

so that 

(92) fs 2m R im ■ ds = [S 2m i? 2m _!)* + 2 m • J S, m _ x R 2m , • A, 



rf*. 



(93) / -S 2m+1 i? 2m+1 -ds — — 2 m+ j [<%„,+! A.^ 4- 2ot . j / ^« -®: 

- £ 

By alternate applications of (92) and (93) any integral of the type 

I S m R m • ds can be reduced to a linear function of expressions of the type 
\_S k -fffc-i]* for values of k descending from m to 1 and of the integral 

~8 S 

(94) fs<,R -ds= Ce'^-ds. 

i. t 

X. Case in which all Values of s are possible. 

For s = — oo and s = + oo (that is, if all values of the residual are 
possible), the integrals in question have simple numerical expressions. 
Namely, whatever integer k may be, 

s* e * = f or * = ± oo ; 
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therefore, 

(95) R k = for * = ± oo ; 

therefore, equations (81) are all satisfied, by (89)— (91). Also, (92) and 
(93) become, for this case, 

+00 -t-00 

I S 2m R 2m • ds = 2 m J S 2m _ ± R im „i • ds, 



+00 +00 



-|-O0 -|-0O 

2 m + 1 J 



by alternate applications of these two formulae, we obtain 

+00 -f-co 

I ^2m ^2m ' « 5 — K , I *->2m-2 -**2m— 2 * "$? 

OO — CO 

-J-co -j-00 

I »2m+l ^2m+l " « 5 ~ 0-.ii I "3m— 1 -"2m-l ' «* > 

—CO — 00 

and again, by repetitions of these, 

+00 -{-00 

J S^R^ ■ ds = _ ' J S R ■ ds, 

00 00 

+00 +00 

J #2m+i -K 2 ,„+i • ds = - '— j S 1 Rr ds. 

CO 00 

Now, by (93), (63), and (68), 

+00 +00 +00 

/ #! #! • ds = j S R ■ ds= fe~* • ds = 4/2V, 

00 — oo 00 

where ir represents the ratio of the circumference of a circle to its diame- 
ter ; therefore, 



+» 



(96) Js^R^-ds^ 



m ! 



(m - £) (m > 

+00 



(97) fs^B^.A^T-^L-fa 



m ! 
(m + iy 



196 



PROCEEDINGS OF THE AMERICAN ACADEMY. 



so that, by (8 


0), 


(98) C, m = 


1 


and 




(99) 


2 „ 



1 (to - i)' 



— AY"» 



2m+l 



V2tt 



for < to < a (excepting w = 1) 

( m + |)M = 

for 1 < to < p. 



It will be observed that 2 „ is less and C 2m+ i is greater than 

for very large values of to we have, by Stirling's formula, 

1 



V2-. 



(100) 



OL 



' \/m 



_ V'2m 



so that C 2 ,„ diminishes and <7 2m+1 increases without limit as to increases 
without limit ; this fact will, perhaps, to a certain extent offset the con- 
clusion that might be drawn from (61) that the values of the P's with 
odd suffixes are likely to be less than those of the P's with nearly the 
same even suffixes (on account of the presence of the odd functions S 
under the integral sign in the former and of the even functions S in the 
latter) in estimating the relative importance of the terms of f(s) as given 
by (78) that contain P's with odd suffixes and P's with even suffixes. 
For the G's with the smaller suffixes we have 



(101) 



C = 



1 



A/2- 



Os — 7= ' 

2 V2*- 



d = 



8 V2n- 

15 



6 
7 = 



5 

16 V^' 
35 



etc., 
etc., 



8 a/2*- 16 a/2 i 

so that, by (27), (78), (68), (51), (63), and (101), 

f(l-ip 4 )(l-2 S 2 +^) 



(102) 



<K*> =£/« = -7== C 1 

P f* a/2 7T 

- I (1 - *P4 + sW (1 - 3 s 2 + s* - -rV* 6 ) + etc. 



i Pa (» - i s 8 ) - f (Ps - A pi) (« - ! * 3 + tV * 6 ) + eto.]e- *• 



XI. Application to ant Particular Case. 

It seems reasonable to suppose that, in any case (whatever the values 

of s and i) f(s) will have small values for values of s beyond the limits 

of possible values, so that the extension of the limits to — co and + » 

will not make much difference in the values of the integrals by which 



(103) {^ 
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the C's are determined ; we may, then, take the values given by (98) 
and (99) as approximate in all cases. 

The restriction of the development (78) olf(s) to a finite number of 
terms is justified only by the fact that the P's that would occur in later 
terms are all negligible, that is that 

f-Pto = ° for a < m, 
t-Psm+i = ° for /S < m, 
to the degree of approximation to which the calculations are restricted. 
In the special case for which s = — ao and s = + oo these equations 
(103) are exactly satisfied by (89)-(91), because the suffixes of the P's 
in (103) are different from those of any of the P's that occur in the 
right members of (79). 

In applying the general method to any particular set of observations 
we have, then, to calculate, by (l)-(7) and (48) (from the observations), 
the values of z , the residuals (£ = z — z ), the p's, the p's, and the P's 
with so many successive suffixes as to satisfy ourselves that the P's with 
suffixes greater than a certain number are all negligible ; the greatest 
even suffix and the greatest odd suffix of the P 's that are not negligible 
are the values of 2 a and 2/3+1, respectively; the values of the P's 
that are not negligible and the approximate values of the corresponding 
C's given by (98) and (99) being substituted in (77) determine an 
approximate form of 9 (s) ; the positive and negative roots of 6 (s) = 
that have the smallest absolute values are approximate values of s and s, 
with which more accurate values of the C's can be calculated by (80) ; 
substituting these values of the C's in (77) we obtain a more exact form 
of 6 (s), from (i) = more exact values of s and s, and from (80) more 
exact values of the C's ; and so we continue until the variations of the 
values of the C's are reduced to negligible quantities ; the final values 
of the C's being substituted in (78) give the most exact form of f(s) 
and, by (27), the most exact form of ib (|) that can be derived from the 
observations. The closeness with which the conditions (81) are satisfied 
by the final values of s and * may be regarded as a test of the goodness 
of the observations, that is of the closeness with which the residuals 
follow the law that they follow most closely. 

It may be observed that there is no case for which a = 1 ; more 
properly, the cases for which a = 1 coincide with the cases for which 

a = 0. 

If P 2 „+i = for 1 = m,f($) as given by (78) is an even function of 
s and the relative frequencies of positive and negative residuals that have 
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the same absolute magnitudes are equal, that is, positive and negative 
residuals are equally probable. 

The calculation of the true values of the G's by (80). (92), and (93) 
will be facilitated by tables of the values of S m R m _y and of the integral 
(94). The values of S m i? m _i will be required for such values of m as 
are suffixes of the P's that are not negligible, which will be small, un- 
doubtedly, in all actual cases ; tables of S„, i?,„-i f° r m = 3, 4, 5, 6 will 
probably be sufficient. For the values of the integral (94) the ordinary 
tables of 

1 ~ 

V- 

can be used ; namely, 



(104) *( r )=— = fe-' 2 -dt 

yir J 



v§ 






/.-*■*- Vi/^-*-Vf.^), 



so that, observing that <J> (— y) = — 4> (y), we have 

a«)/.^-*-»1-[»(^)-^)] ! 

where * and — * are both positive. 

XII. Probable Values. 

We have defined the probable value of the observed quantity to be that 
value that is just as likely to be exceeded as not ; if the probable residual 
£ p and the probable relative residual s p be taken to correspond to this 
probable value z p , £ p and s p will be determined by the equation 
*p s p 

(106) |V (0 ■ *i =//(») - * = * 
) r. 

and z p by 

(107) z p = z + ( p , 

in accordance with (3), so that the arithmetical mean is (in general) not 
the probable value even of a quantity directly observed (as it is assumed 
to be in the method of Gauss). 
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XIII. The Ordinary Case. 
In the ordinary theory it is assumed that 

(108) -P* = for 3 ^ £, 

so that the only P regarded as not negligible is P = 1, which being 
substituted in (77) with the approximate value C — — =, by (101), 
gives, by (51) and (63), 

0<» = c p & = -$=■, 

(s) = has, then, no real finite root, so that 

s = — oo,s = + oo 

and (80) gives again the same value of C as (101), which is, therefore, 
the most exact value possible on the assumption made. We have fur- 
ther, by (78) and (27), 

(109) f( t)= -L = e - 1 j (fi = _J^ e -u* . 

y2ir / i V2 7r 

while (60) gives 

Pk = —j= I s k ■ e * • ds, 

V2 7T J 

that is 

(110) p 2m+ i = for O^ffl 
and 

(111) P2m = 2" (m - !)""> or p 2m = 2» (m - $)<*> • ,«*» 

for < m, so that, far from being negligible for the greater values of 
m, p 2m increases beyond all limit as m increases. 

XIV. Two Other Particular Cases. 

For the purpose of illustrating the variety of forms the general theory 
may take in particular cases, and especially its possible deviation from 

the ordinary theory, we have plotted the curve represented by y = -/(f) 

(the probability-curve) for a = 0, /3 = 1 and for a = 2, ft — 0, using the 
approximate values of C , C 3 , and C± given by (101), for such different 
values of p s and p t (in the two cases, respectively) as to show all the 
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essentially different forms for these pairs of values of a and /3. Figure 1 
shows the locus of the equation 

(U2) y = -i= [1 - i Ps (i - i s 3 )] e _|2 

ft V 2 TT 

for p 3 = 0, 1, 1.648, 2.2, 3, 5, and Figure 2 the locus of the equation 
(113) y = —^= [1 - | (1 - I p 4 ) (1 - 2 s 2 + J S <)] e"* 

jW V 2 TT 

for p 4 = 0, 1, 1.4, 2, 3, 4, 6, 7, 11. The centre of each figure corre- 
sponds to s = and the two extremes to s = ± 4; the horizontal line 
is the axis of abscissae. It will be seen that the ordinates of all the 
curves are practically for values of s whose absolute magnitudes 
exceed 4. The curve for which p 3 — in Figure 1 and that for which 
p 4 = 3 in Figure 2 coincide with the ordinary probability-curve. The 
curve for which p s = 3 in Figure 1 and that for which p 4 = 7 in Figure 2 
touch the axis of abscissae ; these values are, by the footnote on page 173, 
the extreme values of p 3 and p 4 , respectively; each figure contains one 
curve that corresponds to a value beyond these limits, which shows the 
discontinuity of s or s as the curve passes through the form that has 
contact with the axis. While p 4 is essentially positive, as a mean even 
power of the relative residual, p 3 may be negative ; if the sign of the 
value of p 3 that corresponds to either curve in Figure 1 is changed, the 
curve will be reversed, right and left. In the cases represented in 
Figure 1, the maximum relative frequency and the probable value of the 
observed quantity correspond to negative values of s ; in the cases rep- 
resented in Figure 2, the maximum relative frequency corresponds to 
s = for some values of p 4 , while there are two maxima for the smaller 
values of p 4 , the probable value of the observed quantity corresponding 
always to * = 0. 

Because s = 1 corresponds to £ = fi, that is to £ 2 = p 2 , where /> 2 is 
the mean square of the residual, it is evident that at least one of the 
limits s and * must be considerably greater than 1 in absolute magnitude, 
in any case; but the curves plotted show that very little of accuracy is 
gained by considering values of s beyond ± 4 in the cases represented 
by them. 

Clark University, Worcester, Mass. 



